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ABSTRACT 


A  branch  and  exclude  algorithm  for  solution  of  the 

N  N 

"knapsack  problem",  max  \      v,x.  where   J   w,x.    W  and 

1=1  1=1 

x,  =  0,1  ,  is  presented  which  requires  relatively  small 
amounts  of  computer  running  time  and  core  storage  allocation, 
In  addition,  a  branch  and  bound  scheme  is  developed.   The 
branch  and  exclude  method  is  then  compared  to  the  branch  and 
bound  method  and  to  a  branch  and ^batj!d  method  given  by 
Kolesar  [2].   Computational  results  are  given. 
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1.    Introduction 

Stated  in  precise  terms,  the  knapsack  problem  is  one  in 
which  one  wishes  to  select  from  among  a  finite  collection  of 
indivisible  objects  a  subcollection  which  maximizes  a  linear 
function  subject  to  a  linear  inequality  constraint  [2] . 
Dantzig  [4]  gives  the  following  example:   a  person  is  planning 
a  hike  and  has  decided  not  to  carry  more  than  70  pounds  of 
different  items,  such  as  bed  roll,  geiger  counter,  cans  of 
food,  etc.   We  try  to  formulate  this  in  mathematical  terms. 
Let  w.  be  the  weight  of  the  i    object  and  v.  be  its  relative 
value  determined  by  the  hiker  in  comparison  with  the  other 

objects  he  would  like  to  have  on  his  trip.   Let  x.  =  1  mean 

+-  v» 
that  the  i    item  is  selected,  and  x.  =0  mean  that  it  is 

i 

not  selected.   We  express  the  weight  limitation  by 

N 

T   w.x.  <  70,  with  x.  =  0  or  1 ,  and  we  wish  to  choose  the 
.  L ,       l  i  —  l 

i=l 

N 

x.  so  that  the  total  value,   J   v.x.  is  a  maximum. 

1  .  L .  1   1 

i=l 

The  knapsack  problem  has  also  been  used  as  a  model  to 
formulate  and  solve  various  aspects  of  problems  arising  in 
capital  investment  [5],  network  reliability  [6],  capital 
budgeting  [8],  and  in  optimal  methods  of  cutting  stock  [7]. 
Dantzig  also  uses  this  method  in  [4]  to  give  a  delightful 
"proof"  that,  of  all  the  forms  of  marriage  (monogamy,  bigamy, 
polygamy)  monogamy  is  the  best  of  all  possible  relations. 

The  knapsack  problem  may  be  formulated  as  an  integer 
programming  problem  and  solved  by  any  one  of  the  following 
techniques:   (i)  the  cutting  plane  technique,  (ii)  techniques 


employing  parallel  shifts  of  the  objective  function  hyper- 
plane,  (ni)  techniques  based  upon  Boolean  algebra,  and  liv) 
combinatorial  methods  which  enumerate  a  restricted  subset  of 
possible  integer  solutions.   For  detailed  references  on  these 
methods  see  Glover  [1].   An  alternative  formulation  of  the 
problem  may  be  developed  through  the  use  of  dynamic  program- 
ming [9] .   The  technique  to  be  presented  in  this  paper  falls 
into  the  class  of  combinatorial  methods  mentioned  in  (iv) 
above. 

In  combinatorial  methods  of  this  type,  the  computer 
memory  and  running  time  requirements  are  usually  quite  large 
for  problems  involving  many  variables. 

We  present  here  first  a  branch  and  bound  algorithm  that 
is  more  in  the  spirit  of  Land  and  Doig  [3]  than  is  [2] .   We 
then  develop  an  algorithm  which  requires  little  computer 
storage  and  running.   Computational  results  are  then  presented 
for  our  two  algorithms  and  for  Kolesar's  method. 
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2.    Formulation  of  the  Problem 

We  formulate  the  knapsack  problem  as  follows : 

N 
(1)   Maximize    T   v.x. 
1=1 


Subject  to    I      vr.  x.         W, 


N 

w . 
li  — 


i=l 
Where  x.  =  0,  1  for  i=l,2, N 

We  can  assume  without  any  loss  in  generality  that  the 

constants  v.  and  w.  are  integers.   Fractions  can  be  handled 
11^ 

by  multiplying  through  by  a  proper  factor.   Non-positive  con- 
stants may  be  handled  by  the  following  method  which  is  due 
to  Glover  [1] : 

1   if   v.  >  0   and  w.  <  0 

•NT  i  —  i  — 

i)   Let  x.  = 

1      0   if   v.  <  0   and  w,  >  0 
i  —  l  — 

ii)   Next  reformulate  the  problem  as: 

N 
Maximize    >   v.y. 
i=l 

N 
Subject  to    1      w.y.    W  , 
i=l   1  1 

Where   y,  =  0,1  for  i=l,2, N, 

x.   if   v. ,w -    0 

and  where  y .  = 

1-x .   if   v. ,w.  <  0 
l       i   l  — 

The  following  example  should  illustrate  this  technique: 

Maximize     3x,  -  2x2  -  4x3  +  x. 

Subject  to   2x.  -  3x«  +  x_  -  x.  _<  0 


11 


Using  Glover's  procedure  this  problem  is  reformulated 
as  : 

Maximize    3y,  +  2y„  -  1 

Subject  to  2y,  +  3y2    4 
and  the  optimal  solution  if  the  integer  constraint  on  x„  is 
relaxed  is  10/3  units  with  X  =  (1,1/3,0,1). 
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3.    A  Branch  and  Bound  Algorithm 

As  shown  by  Dantzig  [4]  the  optimal  fractional  solution 


to 


N 

(2)   Maximize      T   v.x. 

.  S   11 
i=l 

N 
Subject  to    I      w.x.    W 
i=l 


0  <    x,  <_  1  (i=l,.  .  .  ,N) 


is  given  by 


X.      = 

1 

1      if 

i    <    r 

X.      = 

1 

0      if 

1    >    r 

X        = 

r 

W-J 

i         r 

w . 

i 

w 
r 

where  r  is  the  least  integer  (0    r  :  N)  for  which 

y    w.  >  W  .   If  no  r  exists  we  have  all  x.  =  1.   If 
.  ^  l  —  l 

1  —     r 

x   =0,  we  have  the  optimal  solution  to  (1). 

If  x   is  fractional  the  value  of  the  objective  function 

is  z(l)  =    T   v.  +  v  x   .   We  consider  z(n) ,  the  value  of 
.  £    i    r  r 
i  <  r 

the  objective  function  at  node  n,  as  the  solution  to  (2)  with 
assigned  variables  added  as  constraints.   A  branch  and  bound 
scheme  is  essentially  as  follows: 

1.      Label  node  1  with  z(l).   Go  to  2. 
2(a).   Find  the  terminal  node  with  the  largest  value  of 
z (n) .   This  is  the  node  at  which  the  next  branch- 
ing will  take  place.   Any  node  (^  1)  contains  the 
effect  of  assigning  values  to  variables  and  solv- 
ing (2)  with  the  assigned  values  of  the  variables 
added  as  constraints. 
13 


(b) .   If  the  solution  at  node  n  has  all  integer  var- 
iables, we  have  achieved  an  optimal  solution  to 
(1).   Stop.   If  not,  go  to  3. 
3(a).   Set  n  =  n+1  and  some  unassigned  variable,  say  x,, 
equal  to  0.   Solve  (2)  with  all  assigned  variables 
added  as  constraints.   Label  node  n  with  the  value 
z  (n)  .   Go  to  3 (b)  . 
(b) .   Set  n  =  n+1  and  x,  =  1.   Solve  (2)  with  all 

assigned  variables  added  as  constraints.   Label 
node  n  with  the  value  z(n).   Go  to  2. 
Kolesar's  algorithm  consists  in  taking  as  x ,  ,  in  step 
3(a),  the  unassigned  variable  with  smallest  index  (i.e., 
t  =  i  for  which  v./w.  is  the  maximum  for  unassigned  variables 

x.). 

We  propose  instead  to  follow  [3]  and  to  take  x,  ,  in  step 
3(a),  the  variable  that  is  fractional  at  node  n. 

To  illustrate  our  branch  and  bound  algorithm,  we  solve 

the  problem  given  in  [2] : 

Index 

1 
2 
3 
4 
5 
6 
7 

and  W  =  100.   Solving  (2),  we  obtain  z(l)  =  140  with  x..  =  1, 

x„  =  1,  x3  =  1/2.   We  label  node  1  with  z(l)  =  140.   Since 

x.,  is  fractional  we  branch  from  node  1  and  proceed  to  step 

3(a).   We  solve  (2)  with  x^  =  0.   This  produces  z(2)  =  135 

with  x,  =  1,  x2  =  1,  x3  =  0,  x.  =  1,  x5  =  1/4.   We  label  node 
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w. 

1 

v. 

l 

30 

60 

50 

60 

40 

40 

10 

10 

40 

20 

30 

10 

10 

3 

TTODE  1 


z(1  )=1ij.O 


:tcde  2 


s(2)=l 35 


].v,  ^rj    ^ 


1TO  DE  £ 


(3)=136 


NODE   7     NODE  ij. 


z(6)=133rr 


hod: 


i-o     A 


z(7)=1U 


NODS  5 


zdj.)=120 


z(5)=120 


NODE  9 


z(8)=133 


z(?)=128 


Figure  1#   Complete  tree  for  the  example 
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2  with  z(2)  =  135.   In  step  3(b)  we  set  x~  =  1.   This  produces 
z(3)  =  136  with  x,  =  1 ,  x„  =  3/5,  x_  =  1.   We  go  to  step  2  in 
the  algorithm.   We  see  that  z(3)  is  the  maximum  for  all  ter- 
minal nodes.   Since  the  solution  at  node  3  has  x„  fractional 
we  branch  from  node  3  and  proceed  to  3(a).   We  solve  (2)  with 
x~  =  1  and  x„  =  0  added  as  constraints.   This  produces 
z(4)  =  120,  with  x,  =  1,  x2  =  0,  x3  =  1,  x,  =  1,  x5  =  1/2. 
We  proceed  to  3(b)  and  solve  (2)  with  x~  =  1  and  x~  =  1. 
This  produces  z(5)  =  120  with  x.  =  1/3,  x2  =  1 ,  x,  =  1.   Re- 
turning to  step  2  we  see  that  z(2)  is  the  maximum  for  all  ter- 
minal nodes.   Since  the  solution  at  node  (2)  has  x,.  fractional 
we  must  branch  from  node  2.   The  method  continues  easily  and 
the  complete  tree  is  shown  in  Figure  1.   The  optimal  solution 
is  given  at  node  8  with  z(8)  =  133,  x.  =  1 ,  x„  =  1,  x.  =  1, 

x_  =  1  and  all  other  x.  =  0. 
7  i 

The  algorithm  achieved  solution  with  the  generation  of 
9  nodes.   The  same  problem  was  solved  in  [2]  and  required  15 
nodes.   Further  computational  experience  is  presented  in 
Section  5. 
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4.    A  Branch  and  Exclude  Algorithm 

Branch  and  bound  techniques  normally  require  large 
amounts  of  computer  storage  because  of  the  necessity  to  store 
information  for  all  terminal  nodes.   In  this  section,  we  pre- 
sent a  branch  and  exclude  algorithm  for  solving  (1)  that  re- 
quires little  computer  storage. 

Our  branch  and  exclude  algorithm  first  finds  an  obvious 
integer  solution  to  the  constraints  of  (1) .   This  solution 
is  a  lower  bound  to  the  optimal  solution.   We  develop  a  branch 
of  a  tree  and  explore  each  part  of  the  branch  until  the  lower 
bound  is  reached  or  until  a  new  feasible  solution  is  found 
that  represents  a  larger  lower  bound.   We  then  back-track 
and  develop  new  branches  of  the  tree  developing  possibly 
larger  lower  bounds.   Further  branching  is  excluded  when  the 
lower-bound  is  reached.   The  algorithm  stops  when  all  new 
branches  are  excluded.   The  only  information  that  is  stored 
is  the  current  lower  bound  solution  and  the  branch  routing. 
At  the  end,  the  lower  bound  solution  is  optimal. 

Define  [x]  as  the  greatest  integer  less  than  or  equal 
to  x. 

Define  S (xfl ,x, , . . . ,x-J  as  the  current  lower  bound  solu- 

N 
tion,  where  the  x.  are  all  given  and  x~  =   )   v.x.. 

1=1 
Define  X (x. ,x„ , . . . ,x^)  to  indicate  assigned  variables. 

A  value  x.  =2  (or  any  number  not  equal  to  zero  or  one)  in- 
dicates the  variable  is  unassigned.   An  assigned  variable 
will  have  the  value  zero  or  one. 

Define  R(L)  as  the  index  of  the  L  assigned  variable. 
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The  solution  to  (2) ,  with  L  assigned  components  of  X 

added  as  constraints,  is 

x.  =  1   if   i  <  r,  i  j-   R(j)   (j  =  1,...,L) 

x.  =  0   if   i  >  r,  i  i-    R(j)   (j  =  1,...,L) 

(3)      x  ,  ,,  =  0  or  1  (j  =  1,...,L)  depending  on  the 
R  \  J  i 

assignment 

L 
x  =  (W  -  l     w    xR  .   -  I        w  )/w 

1=1  ieM(L) 

Z(L+1>  =   J   VR  (i)KR      +   X     V   ♦  V  X 

1=1  ieM(L) 

where  the  set  M  is  given  by 

M(L)  =  (i|i  <  r,  i  ±    R(j)  (j  =  1,...,L)} 

and  r  is  the  least  integer  (0  _<  r   N)  for  which 

L  L 

I  w.  +  w   >  W  -  I      w  mx  m   . 

ieM(L)   X     r        i=l   RU)  RU) 

If  no  r  exists  we  have  all  x.  =  1  for  i  ^  R(j)  (j  =  1,...,L). 

A  lower  bound  to  the  solution  of  (1)  is  given  by 

L 

i=l   R(l)  R(l)    ieM(L)   X 

The  algorithm  follows: 

1.  Set  L  =  1  and  all  components  of  X  to  two.   Go  to  2. 

2.  Solve  (2) .   If  'the  solution  is  all  integer  we  have  the 

optimal  solution.   Stop.   If  not  the  solution  is  x,  =1, 

x0  =  1.....X   ,  =1  and  x   is  fractional.   We  calculate 
2     '    '  r-1  r 

x„  =   J   v.  and  form  S (xfl ,1,1, . . . ,1 ,0 ,0 ,0)  as  a  lower 
i<r   1 

bound  to  the  optimal  solution,  where  the  zero  components 
in  S  represent  x.  =  0,  i    r.   Set  R(l)  =  r  and  the  r 
component  of  X  to  zero.   Go  to  3. 
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3(a).   Solve  (2)  with  the  L  assigned  components  of  X  added 
as  constraints. 

We  obtain  (3).   If  [Z(L  +  1)]  <  xQ  to  to  4.   If 
Z(L  +  1)  >  x0  and  we  have  an  integer  solution,  take 
xQ  =  Z(L  +  1) ,  form  a  new  S  (xQ  ,xx  ,x2  ,  .  .  .  ,xN)  from  (3),  and 
go  to  4.   If  [Z (L  +  1) ]  >  xQ  and  we  do  not  have  an  integer 
solution,  go  to  3(b). 

(b) .   If  x(L)  >  xQ ,  take  xQ  =  x(L)  and  form  a  new 
S(xQ,x1,. . . ,xN)  from  (3)  with  xr  =  0 .   In  any  case,  Set 
L  =  L  +  1,  take  R(L)  =  r,  and  set  the  R(L)  component  of  X  to 
zero.   Go  to  3 (a) . 

4(a).   If  the  R(L)  component  of  X  is  equal  to  zero  change  the 
component  to  one  and  go  to  3(a).   If  the  R(L)  component  of 
X  is  one  go  to  4(b). 

(b) .   If  L  =  1  the  optimal  solution  is  S (xQ ,x, , , . . ,xN) . 
Stop.   If  L  /  1  change  the  R(L)  component  of  X  to  two,  set 
L  =  L-l ,  and  go  to  4(a). 

This  completes  the  algorithm,  z (L)  represents  the  value 
of  the  objective  function  at  the  L    level  of  a  branch.   R(L) 
represents  a  routing  of  assignments  along  the  branch.   Only 
one  branch  is  studied  at  a  time  with  preference  given  to 
assigning  the  value  of  zero  to  the  fractional  variable.   This 
allows  lower-bounds  to  be  achieved  more  rapidly.   The  only 
permanent  storage  information  required  is  the  current  lower 
bound  solution  S,  the  assigned  variables  R(L)  at  level  L  of 
the  branch,  and  the  assignment  vector  X. 
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In  the  problem  given  above  the  algorithm  as  viewed  in 
Figure  1,  goes  as  follows: 

1.  At  node  one  the  optimal  solution  to  (2)  is 
fractional.  x~  =  120  . 

2.  Branch  to  node  two,  assigning  x^  =  0. 

3.  Branch  to  node  six,  assigning  x_  =  0. 

4.  Branch  to  node  eight,  assigning  x-  =  0 .   This 
produces  an  all  integer  result  with  x~  =  133, 

X.  =  X  ,  Ky    =  I  ,  X  *  =  I  ,  X—  =  X . 

5.  Back-track  to  node  six;  branch  to  node  9, 
assigning  xg  =  1. 

6.  Back-track  to  node  6,  removing  the  assign- 
ment on  x,. 

6 

7.  Back-track  to  node  2;  branch  to  node  7, 
assigning  X-  =  1. 

8.  Back-track  to  node  2,  removing  the  assign- 
ment on  x_. 

9.  Back-track  to  node  1;  branch  to  node  3,  assign- 
ing x3  =  1.   Calculate  x(l)  =  100;  do  not 
change  the  current  lower  bound. 

10.  Branch  to  node  4,  assigning  x?  =  0. 

11.  Back-track  to  node  3;  branch  to  node  5,  assign- 
ing x2  =  1 . 

12.  Back-track  to  node  3,  removing  the  assignment 
on  x„ . 

13.  Back-track  to  node  1;  end. 

14.  The  optimal  solution  is  given  at  node  8. 

In  this  example,  the  total  tree  investigated  has  the  same 

number  of  nodes. 
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5.    Computational  Results 

We  have  programmed  the  two  branch  and  bound  schemes  and 

the  branch  and  exclude  algorithm  in  Fortran  IV  and  have  run 

test  problems  on  the  IBM  360/67.   A  program  for  the  branch 

and  exclude  method  is  given  in  Appendix  A. 

The  coefficients  v.  and  w.  were  generated  as  random 

1      1      ^ 

integers  and  the  three  methods  were  then  tried  on  the  same 

problem.   If 

BK  =  Kolesar's  branch  and  bound  algorithm, 

BB  =  the  branch  and  bound  algorithm  presented  here, 

BE  =  the  branch  and  exclude  algorithm  presented  here, 

we  obtained  the  following  results  for  the  number  of  nodes 

generated  and  the  time  to  achieve  solution: 

TABLE  1.   Comparison  of  Methods 


w 

N 

20 

30 

40 

50 

20 

30 

40 

50 

BE 

28.8 

28.0 

28.8 

35.7 

61 

87 

120 

164 

25 

BB 

26.4 

20.6 

26.6 

28.6 

201 

282 

395 

494 

BK 

56.2 

49.6 

48.0 

50.0 

259 

345 

447 

552 

BE 

34.5 

38.7 

36.1 

54.5 

75 

111 

145 

235 

50 

BB 

38.6 

48.2 

.  42.0 

51.8 

236 

361 

453 

581 

BK 

62.2 

78.4 

78.0 

92.6 

276 

434 

553 

723 

BE 

19.8 

39.0 

27.8 

18.5 

58 

116 

128 

135 

75 

BB 

31.0 

30.2 

27.2 

25.0 

222 

312 

404 

485 

BK 

62.6 

101.2 

84.4 

61.6 

280 

523 

578 

616 

BE 

7.2 

24.3 

46.8 

29.0 

36 

97 

181 

163 

100 

BB 

16.6 

24.8 

64.0 

50.4 

191 

301 

560 

598 

BK 

[56.6 

65.4 

136.6 

134.2 

273 

403 

806 

978 

Average  number  of  nodes  generated.   Average  time  in  Milli- 
seconds.  (10  problems  run  to  obtain  each  average) 
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In  these  test  problems  superiority  of  BB  and  BE  over  BK 
is  clearly  apparent  with  regard  to  number  of  nodes  generated. 
As  for  time  considerations, the  branch  and  exclude^  method 
shows  to  be  markedly  faster  over  the  other  two  methods,  while 
BB  is  faster  than  BK.   Various  other  problems  than  those 
listed  above  were  run,  and  in  general,  the  BE  method  achieved 
solution  in  about  one-half  to  one-third  the  time  required  by 
the  BK  and  BB  methods. 

The  branch  and  exclude  method  was  tried  on  several  prob- 
lems with  5000  variables  and  solution  was  achieved  in  all 
cases  in  approximately  4  minutes.   The  other  two  methods 
failed  because  of  storage  limitations.   In  general,  the  branch 
and  exclude  algorithm  presented  here  can  rapidly  solve  ex- 
tremely large  problems. 


22 


BIBLIOGRAPHY 

1.  Glover,  F. ,  "A  Multiphase  Dual  Algorithm  for  the  Zero-one 

Integer  Programming  Problem",  Operations  Research,  Vol. 
13  (1965) ,  pp.  879-919. 

2.  Kolesar,  P.  J.,  "A  Branch  and  Bound  Algorithm  for  the 

Knapsack  Problem",  Management  Science,  Vol,  13  (1967), 
pp.  723-735. 

3.  Land,  A.  H.,  and  Doig,  A.  G. ,  "An  Automatic  Method  of 

Solving  Discrete  Programming  Problems",  Econometrica, 
Vol.  28  (1960)  ,  pp.  497-520. 

4.  Dantzig,  G.  B.,  "Discrete  Variable  Extremum  Problems", 

Operations  Research,  Vol.  5  (1957)  ,  pp.  256-277. 

5.  Hansmann,  F. ,  "Operations  Research  in  the  National  Plan- 

ning of  Underdeveloped  Countries",  Operations  Research, 
Vol.  9  (1961)  ,  pp.  203-248. 

6.  Kolesar,  P.,  "Assignment  of  Optimal  Redundancy  in 

Systems  Subject  to  Failure",  Columbia  University, 
Operations  Research  Group,  Technical  Report,  1966. 

7.  Gilmore,  P.  C. ,  and  Gomory ,  R.  E.,  "A  Linear  Programming 

Approach  to  the  Cutting  Stock  Problem-Part  II",  Opera- 
tions Research,  Vol.  11  (1963) ,  pp.  863-888. 

8.  Cord,  J.,  "A  Method  for  Allocating  Funds  to  Investment 

Projects  when  Returns  are  Subject  to  Undertainty" , 
Management  Science,  Vol.  10,  No.  2  (1964),  pp.  335-341. 

9.  Bellman,  R.  E. ,  and  Dreyfus,  S.  E.,  Applied  Dynamic 

Programming,  Princeton  University  Press,  1962, 
pp.  27-31. 


23 


APPENDIX  I 

FORTRAN  IV  Computer  Program  for  the  Brancv  and  Exclude 

Method 

LEVEL  It  MOO  0  MAIN  OATE  *  68051 

IMPLICIT  REALCY-Z, S ), INTEGER <A-X > 

INTE6Eft*2  V,W,OROER,X,SVECTiR,F,G,H 

COMMON  Vt 100) tW< 100 »fN*X< 100  1 tOROER< 100) tWNOT, 
*Y  VEC < 100 ) , SFLAG, XNOT , SVECT< 100 ) ,L,R  < 100  > , F ( 100 ) , 
*G<100) 

DIMENSION  H<100> 
C 

C  V,W,R,X  AND  L  ARE  AS  STATED  IN  THE  TEXTUAL  EXPLANATION  OF 
C     THIS  ALGORITHM 

C  WNOT  IS  lH€r   WEIGHT  CONSTRAINT,  W 
C  N  IS  THE  NUMBER  OF  VARIABLES 

C  SVECT  IS  THE  SOLUTION  VECTOR,  S,  EXCEPT  THAT  X-SUB-ZERO 
C     IS  KEPT  AS  XNOT 

C  SFLAG  IS  THE  NUMBER  ASSIGNED  TO  A  PARTICULAR  NOOE 
C 

BEAD  101, WNOT, N 

101  FCRMAT<2I5) 
DO  102  1  =  1, N 

READ  ioi,vm,*<  n 

102  CONTINUE 

DO  103  1  =  1, N 
R(  I)=0 
X(I)=2 
SVECTU)=0 

103  CONTINUE 
*NOT=Q 
SFLAG=0 
L*l 

C 

C  THE  VARIABLES  ARE  HERE  PLACED  IN  ORDER  OF  DESCENDING 

C     VSUBI/WSUBI 

C 

CALL  RANK 
C 

C  THE  CPTIMAL  SOLUTION  IS  NOW  FOUND  BY  THE  SOLVIT  ROUTINE 
C 

CALL  SCLVIT 
C 

C  AT  THIS  POINT  THE  VARIABLES  ARE  REARRANGED  INTO  THEIR 
C     ORIGINAL  ORDER*  AND  THE  RESULTS  PRINTED 
C 

DO  91  T-UN 

FCCRDERU)  )  =  SVECT(  I) 

G(CRDERM)  )=V(  I) 

H(CRDER( I) )=W( I ) 

91  CONTINUE 
PRINT  70, XNOT 

70     F0PMAT<1H1, 20X, 'SOLUTION  =  •, 15, •  UNITS1) 

PRINT  90 
90     FORMAT* 20X  r»  I '  rV9X,  *X<  I )  Wi^rX,  »  V<  I )  •  ,  H.X,  •  W(  I  M  ) 

CO  92  1=1, N 

PRINT  93,I,F(I),G( I),H(I) 
93     F0RMAT(3X,4(17X,I3)) 

92  CONTINUE 
STOP 
END 


-2± 


IMPLICIT  RC  ALiV-Zy^hr  INTEGER  <A-*» 

INTEGER*2  V, W, ORDER, X,  SV EC TfR,F,G»H 
COMMON  V ( 100 ) . W ( 100 ) , N f X ! 100 ) , OROER ( 100 ) , WNOT « 
*YVEC<1001,SFLAG,XNOT,SVECT<100)»L.R<1001,F<100). 
*G< 100) 
C 

C  OROER  IS  THE  VECTOR  WHICH  WILL  CONTAIN  THE  OROER  IN  WHICH 
C     THE  VARIABLES  ARE  REARRANGED 
C 


DO    1    1*1, N 
ORDER (I)*! 

z»v(ii 

YVECU)»Z/Y 

1  CONTINUE 
C 

C    THE    RAN*   OROER    IS    NOW    DETERMINES   ANO    STORED    IN    OROER 

K*0 

2  K*K*l 
J*K*1 

ZAP*YVECtK) 
M*K 

DO    3    I=J,N 

*F<VVE€44KLE.ZAP)GG    TO    3 
ZAP*YVEC(I> 
H»I 

3  CONTINUE 
YVEC<M)=YVECIK) 
YVEC(K)=ZAP 
PIP*OROER(M» 
CRDER(M)=ORDERIK) 
CRDER<KI=PIP 
IFU.NE.MGO    TO    2 

C 

C  THE  VARIABLES  ARE  NOW  REARRANGED 


C 


DO  4  1*1. N 
FMI-VCI) 

G(I)=W(I) 

vun i i nuc 

00  5  I*itN 

vm  =  F<ORDERCI)  ) 

*( I )=G<0RDER<  I)  ) 

CONTINUE 

RETURN 

END 


SUBROUTINE  SOLVIT 


IMPLICIT  RCAL(Y-Zt^H+NTEGER<A-X ♦ 
INTEGER*2  V , W ,  ORDER , X t SVECT, R , F ,G,H 
COMMON  V ( 100 ) , W < 100 ) , N t  X ( 100  I .OROER ( 100 1 , WNOT • 
*YVEC(100».SFLAG.XN0T,SVECTI100).L.R(100»,F(100). 
♦  GUOO) 
C 

C  SUM  IS  THE  VALUE  OF  THE  SUM  OF  WSUBI*XSUBI  FOR  ALL  I 
C  TOTALIS  THE  VALUE  OF  THE  SUM  OF  VSUBI*XSUBI  FOR  ALL  I 
C  VVEC  IS  A  VCCTQR  WHICH  HQIOS  THE  VALUE  OF  THE  VARIABLES 
C     AS  IACM  WOE  IS  GENERATED 
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10       Wetf.W'" 


11     CONTINUE 
3      SUM-0 

TOTAL*0 

SFLAG-Sf4.AG*i 
C 

C  THE  VARIABLES  WHICH  HAVE  BEEN  CONSTRAINED  TO  BE  ONE  ARE 
C     HERE  CONSIOEREO 
C 

00  12  I*i,N 

IF<X<  D.NE.DGO  TO  12 

YVECU)«1 

sum»sum»wu> 
tctal-total  +  vu) 
if<sum-wn0t>12,14,20 
12        continue 

C 

C  IF  THE  CONSTRAINT  EQUATION  HAS  NOT  BEEN  VIOLATEO, 

C     UNASSIGNEO  VARIABLES  ARE  NOW  CONSIOEREO 

C 

00  13  I*ltN 

IF(X(  D.LE.DGO  TO  13 

YVEC(1I=1 
C 

C  THE  INITIAL  LOWER  BOUND  VECTOR, S,  IS  BEING  FORMED  BY  THE 
C     NEXT  CARD 
C 

IF(SFLAG,EQ.l)SVECTm  =  l 

SUM«SU*HMI» 

TCTAL=TOTAL  +  VUI 

IF(SUM-WN0TI13. 14,15 

13  CONTINUE 
C 

C  IF  STATEMENT  14  IS  REACHED,  AN  INTEGER  SOLUTION  HAS  BEEN 

C     FCUND 

C 

14  CAtt~STOR€(  TOTAL* 
IFISFLAG.EO.DRETURN 

C 

C  THE  VALUE  ASSIGNED  TO  THE  FRACTIONAL  VARIABLE  IS  NOW 

C     CHECKEO,  ANO  BRANCHING  INSTRUCTIONS  GENERATED 

16  C=X<R<L>) 
IF<Q.EQ.O)GO  TO  17 

20     XCRCLH-2 
L=L-l 

IFU.EQ.DRETURN 
GO  TO  16 

17  X(R<LM  =  1 
GO  TO  10 

C 

C    IF    STATEMENT    15    IS    REACHEO,    A    NON-INTEGER    SOLUTION    HAS 

Co  c  i*m     rni  ihin 
Been    ruuNU  _ 

C 

15  SUM*SUM-WU) 
TOTAL*TOTAL-V(I) 
YVECUI-0 

C 

C  THE  NEXT  CARD  CHECKS  TO  DETERMINE  IF  A  NEW  LOWER  BOUND 

C     CAN  BE  ENTEREO 

mMRfLH«CQ.l.ANP. TOTAL. GT.XNOTKALL  STORE«TOTAL> 

C 

C  FORMULATION  OF  THE  INITIAL  LONER  BOUNO  IS  COMPLETEO  BY 

C     TME  NEXT  THREE  CARDS 

C 
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IF4fFLAG.NE.llG0  TO  79 
XNOT*TOTAL 


C 

C  A  DETERMINATION  IS  NOW  MAOE  AS  TO  WHETHER  THE  GREATEST 

C  INTEGER  CONTAINEO  IN  THE  VALUE  OF  THE  PRESENT 

C  NGN-INTEGER  SOLUTION  IS  EQUAL  TO  THE  VALUE  OF  THE 

C  LOWER  BOUNO  SOLUTION 

C 

79     Z*WNOT-SUM 

v«wm 

¥V€C-4ff»Z/Y 

IFCSFLAG.NE.  UGQ    TO    2 

ZSOL*TOTAL*Vm*YVEC<  I) 

SINT=TQTAL+(WNOT-SU*)*Vm/WU) 

IF4  5fNT.EQ.XNOT  I RETURN 

GO    TO    6 
2  Z»TOTAL*Vm*YVECUI 

SINT*TOTAL+< WNOT-SUMI*V( I >/W< I) 
6      IFtSINT.L-€^XNOTlGQ  TQ    1^ 
C 

C  IF  THIS  POINT  HAS  BEEN  REACHEOt  THE  LARGEST  INTEGER 
C     CONTAINEO  IN  THE  VALUE  OF  THE  CURRENT  NON-INTEGER 
C     SOLUTION  IS  LARGER  THAN  THE  VALUE  OF  THE  LOWER  BOUNO 
C     SOLUTION.  AND  SO  THE  FRACTIONAL  VARIABLE  IS  ASSIGNEO 
C     A  ZERO  VALUE 
C 


IFCYVECUI.LT. 1. AND. YVECUI.GT.OIGO    TO    19 

18  CONTINUE 

19  L*L*1 

5      CONTINUE 
ML)*  J 
X(J)«0 
GO  TO  10 
END -    - 


Implicit  REAL(v-l>>)f Integer i a -hi 

INTEGER*2  V 9W, ORDER. X, SVECT.RtF.G,H 


COMMON  V(  100  >•  WMOO )  t  N.X(  1001  .ORDER  <  100  >,WNOT. 

*YVECTlOO),SFLAG,XNOTfSVECTC100»,L,R(100»,F(100l. 

♦G^lOO) 
C 

C  IF  THIS  ROUTINE  HAS  BEEN  CALLEDt  TOTAL  IS  THE  VALUE  OF  AN 
C     INTEGER  SOLUTION 

C TOTAL  WILL  BE  COMPARCD  TO  XNOT.  1HE    VALUE  OF  THE 

C     CURRENT  LOWER  BOUND  SOLUTION.  AND  THE  SMALLER  OF  THE 
C     TWO  SOLUTIONS  WILL  BE  REPLACEO  BY  THE  LARGER 


C 


IF(TGTA4..LE^XN0TIRETURN 
XNCT-TOTAL 
00  4  1*1. N 
SVECT(U*YVEC(II 
CONTINUE 


RETURN 
END 
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APPENDIX  II 


Logic  Diagram  for  the  SOLVIT  Subroutine  Contained 

in  Appendix  I 


Enter  SOLVIT  Subroutine 


Find  and  store  the  largest  lower 

bound  solution,  S(xA,x, ,...x  ) 

0   1      n    I 


I 


Can  an  optimal  (not  necessarily 
a   integer) solution  be  found  which 
satisfies  all  constraints? 


Yes 


No 


Is  the  optimal  solution  an 
integer  solution? 


No 


Yes 


Is  z(L  +  1)    xQ? 


Yes 


No 


Is  [z(L  +  L) ]  <  xQ? 


No 


....  Yes 


Set  R(L)  =  0 

Set  L  =  L  +  1 


Replace  S(xn,x, ,...x  ) 

with  the  current  integer 
solution 


Is  R(L)  =  1? 


Set  R(L)  =  1 


Yes 


Is  L  =  1? 


Yes 


Set  R(L)  =  2 
Set  L  =  L  -  1 


Return 
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